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II.  ABSTRACT 

It  is  proposed  that  the  usual  energy  deposition,  e.g. ,  electron  beam,  experiments,  as 
well  as  the  usual  shock  introduction ,  e.g. ,  flyer  plate,  experiments  do  not  yield  a  true 
Gruneisec  parameter  of  porous  materials.  It  is  proposed  that  in  a  porous  material  stress 
relief  occurs  wltum  the  time  that  energy  is  being  deposited,  or  before  the  final  state  is 
readied  in  a  flyer  plate  experiment.  Both  types  of  experiments  are  analyzed  and  correlated 
and  agree  with  the  ideas  put  forwar  '  here.  Energy  deposition  experiments  of  Shea, 

Mazella,  and  Avrami  (porous  PETN),  and  fiyer  plate  experiments  of  Boade  (sintered 
porous  copper) ,  are  considered. 

Wave  propagation  in  porous  materials  is  considered  in  some  detail,  and  it  is  concluded 
that  in  such  materials  macroscopic  measurements  can  not  always  be  used  to  determine 
microscopic  parameters.  It  is  proposed  that  a  porous  material  acts  as  a  diffraction 
grating  with  respect  to  a  shock  wave,  and  that  unidirectional  one-dimensional  strain  does 
not  hold. 


1473 


UNCLASSIFIED 
Security  Claaaifi  cation 


1 


1  ORIGINATING  ACTIVITY:  Enter  the  a«ne  end  address 
?f  th*  contractor,  subcontractor,  grantee.  Of)  art  pent  of  De¬ 
fense  activity  or  other  organization  ( corporate  author)  issuing 
»h*  report. 

2a.  REPORT  SFCUWTY  CLASSIFICATION:  Enter  the  over- 
alt  security  classification  of  the  repost.  Indicate  whether 
••Restricted  I rata**  is  included.  Harking  is  to  be  in  accord-  , 
ante  *nth  appropriate  security  regulations. 

2b.  GROUP:  Automatic  downgrading  Is  apeci'led  in  DoD  DP  | 
rective  3200. 10  and  Artred  Forces  Industrial  Manual.  Enter 
the.  group  number.  Also,  when  applicable,  show  that  optional  j 
markings  have  been  used  for  Group  3  sad  Group  4  as  author-  i 
»zed. 

3.  REPORT  TITLE:  Enter  the  complete  report  title  in  all 
capital  latter e.  Titles  in  all  cases  should  be  unclassified. 

If  t  meaningful  title  cannot  be  aelectad  without  cUeeifica* 
non.  show  till#  classification  in  all  capUaia  in  parenthesis 
immedisvely  following  the  title. 

4.  DESCRIPTIVE  NOTES:  If  appropriate,  enter  the  type  of 
teport.  e.g,,  Interim,  progress,  summary,  annual,  or  final. 

Give  the  inclusive  dates  when  a  specific  reporting  period  is 
covered. 

5.  AUTHORS).  Enter  the  namefo)  of  sutboKs)  as  shown  on 
or  tn  the  report.  Enter  last  name,  first  name,  middle  initial. 

If  military,  «ho  r  ra.ik  and  branch  of  service*  The  name  of 
the  principal  author  ia  an  abaolute  minimum  requirement. 

6.  REPORT  DATE  Enter  the  date  of  the  report  as  day, 
month,  vear,  or  month,  year.  If  more  than  one  date  appears 
on  the  report,  use  date  of  publication. 

7*  TOTAL  NUMBER  OF  PAGES-  The  total  page  count 
should  follow  normal  pagination  procedures,  i^e.,  enter  the 
number  of  p*«gec  containing  iruo*mat:on. 

7£>  NUMBER  OF  REFERENCES  Enter  the  total  number  of 
referen  es  cited  in  the  report. 

Sa  CONTRACT  OR  GRANT  NUMBER:  If  appropr.ate,  enter 
the  «pp1ic*h'e  number  of  the  contract  or  grant  under  which 
the  report  wa«  written. 

8.’;,  &  ,  U  8J  PROJECT  NUMBER  Enter  the  appropriate 
military  department  idenMf. cation,  such  os  project  number, 
nubprcject  number,  system  numbers,  task  number,  etc. 

Qa  ORIGIN  A  ICR’S  REPORT  NUMRER(S).  Enter  the  offi¬ 
cial  report  number  by  which  the  document  will  be  identified 
and  controlled  the  originating  activity.  This  number  must 
I  be  unique  to  this  report. 

Qb  OTHER  REPORT  NUMBER* S)  If  the  report  has  been 
i  aaMgnjd  any  ot'.er  report  rumbe.s  (either  by  the  originator 
o'  bv  fhe  tpon^or;,  alto  enter  thii  number* a). 


10.  AVAIL  ABILIT //LIMITATION  NOTICES  Ewer  any  lim 
ita'iocs  on  furthay  d,»aeminaUo*f  of  tha  report.  otirnr  thau  the 
u^oh^  by  tKority  cl.. aificatioti,  „in|  rtwd.it1  ,l.tw»nl 
•uch  art 

(.)  "Qualified  riqwrtKi  tsay  obtain  copies  of  thi* 
report  born  DDC.” 

(2)  "Foreign  t two unr mint  amt  disacmination  of  thia 
retort  by  DDC  1*  not  euthoriraa.  *' 

(3)  "U.  S.  Govrttnwrt  agencies  may  obtain  copiaa  of 
tlua  report  directly  from  DDC.  Othae  qualified  DDC 
aanre  .hall  request  through 

*» 

(4)  "U.  S.  military  aganciaa  may  obtain  copiaa  of  tbit: 
report  dtre<-tly  bom  DDC  Other  qualified  uaeta 
ahall  requeal  through 


(5)  "All  dimtribution  of  thia  report  is  controlled.  Qu»l- 
ifted  DDC  uaare  ahall  raqueat  through 


If  the  report  haa  been  iumiabed  to  the  Office  of  Technical 
Service..  Department  of  Commerce,  for  eale  to  the  public.  Indi¬ 
cate  this  fact  and  enter  the  price,  tf  known. 

1L  SUPPLEMENTARY  NOTES:  Uee  for  additional  enplana- 
tory  notea- 

12.  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  departmental  project  office  or  laboratory  sponsoring  'pay¬ 
ing  lor)  the  reaearch  and  development  Include  addreas. 

13  ABSTRACT  F.nter  an  abstract  giving  a  brief  and  factual 
summery  of  the  document  indicative  of  the  report,  even  'hough  | 
it  may  also  appear  elsewhere  in  the  body  of  the  technics-.  ; 
port  If  edditionet  space  is  required,  e  continuation  sheet  , 
shfcll  be  attached. 

It  IS  highly  desirable  that  the  abstract  of  classified  re¬ 
ports  be  unclassified.  Each  paragraph  of  the  abstract  shell 
end  with  an  indication  of  the  military  security  cUssilhcation 
of  tne  informauon  in  the  paragraph,  represented  ns  (TS),  (i), 

(C)  or  (U) 

There  m  no  limitation  on  the  length  of  the  abstract.  How¬ 
ever  the  suggested  length  is  from  150  to  .25  wordt- 
14,  KEY  WORDS  Key  words  are  technically  meaningful  terms 
or  short  phrases  that  characterize  a  report  end  may  be  used  as 
index  entries  for  cataloging  the  report  Key  words  must  be 
selected  so  that  no  security  classification  is  required  Iden- 
'iers,  such  as  equipment  model  designation,  trade  name,  mili¬ 
tary  project  code  name,  geographic  location,  may  be  used  as 
Wev  words  but  will  be  followed  bv  an  indication  of  ’echnical 
context  The  assignment  of  links,  rules,  and  weights  is 
opt.onal  _ 


“a5S'  jgji4,ywt>H?  gyyu«,-^5*?s^* 


PICATINNY  ARSENAL  TECHNICAL  REPORT  |255 
•* 

THE  GRUNEEEN  CONSTANT  OF  POROUS  MATERIALS  IN  ENERGY 

DEPOSITION  MATERIALS 

« 

BY 

PAUL  HARRIS 
AUGUST  1971 


Approved  for  public  release-  distributor!  unlimited. 


CONCEPTS  AND  EFFECTIVENESS  DIVISION,  NED 
PICATINNY  ARSENAL 
DOVER,  NEW  JERSEY 


1 


ACKNOWLEDGMENT 


The  author  wishes  to  thank  Mr.  Louis  Avrami  of  Picatinny  Arsenal's 
Explosive  Laboratory  for  suggesting  the  problem  to  him,  and  for  many 
helpful  discussions. 


TABLE  OF  CONTENTS 


PAGE  NO. 


List  of  Illustrations  4 

List  of  Tables  5 

Abstract  .  g 

Summary  7 

Section  L  Introduction  g 

IL  The  Griineisen  Constant  10 

UL  Microscopic  Effects  and  the  Griineisen 

Constant  13 

IV.  Shock  Propagation  in  Porous  Material  s  „  19 

V.  Microscopic  Effects  on  Shock  Propagation  20 

References  23 

Distribution  List  25 


3 


LIST  OF  ILLUSTRATIONS 


FIGURE  NO.  PAGE  NO. 

1  An  Averaging  Model  for  Effective  Pressure  Due 

to  Energy  Deposition  in  a  Porous  Solid  15 

2  Shock  Velocity  versus  Particle  Velocity  for 

Sintered  Porous  Copper  17 

3  Model  for  a  Porous  Solid  19 


4 


LIST  OF  TABLES 


TABLE  NO.  PAGE  NO. 


I 

Electron  Mean  Free  Paths  in  Some  Materials 
as  a  Function  of  Energy 

11 

n 

Energy  Deposition  Properties  of  PETN 

12 

m 

Grttneisen  parameters  for  Sintered  Porous 

Copper 

18 

5 


ABSTRACT 


It  is  proposed  that  the  usual  energy  deposition,  e.g. ,  electron  beam,  experi¬ 
ments,  as  well  as  the  usual  shock  introduction,  e.g.,  flyer  plate,  experiments  do  not 
yield  a  true  Gruneisen  parameter  of  porous  materials.  It  is  proposed  that  in  a  porous 
material  stress  relief  occurs  within  the  time  that  energy  is  being  deposited,  or  be¬ 
fore  the  final  state  is  reached  in  a  flyer  plate  experiment.  Both  types  of  experiments 
are  analyzed  and  correlated  and  agree  with  the  ideas  put  forward  here.  Energy  deposi¬ 
tion  experiments  of  Shea,  Mazella,  and  Avrami  (porous  PETN),  and  flyer  plate  experi¬ 
ments  of  Boade  (sintered  porous  copper),  are  considered. 

Wave  propagation  in  porous  materials  is  considered  in  some  detail,  and  it  is 
concluded  that  in  such  materials  macroscopic  measurements  can  not  always  be  used  to 
determine  microscopic  parameters.  It  is  proposed  that  a  porous  material  acts  as  a 
diffraction  grating  wiih  respect  to  a  shock  wave,  and  that  unidirectional  one-dimensional 
strain  does  not  hold. 
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SUMMAHY 


In  this  report,  we  have  looked  at  some  of  the  details  associated  with  shock 
propagation  in  porous  materials.  We  have  inquired  into  effects  resulting  from 
shock  generation  by  "instantaneous"  energy  deposition  or  by  shock  introduction, 
e.  g. ,  via  a  flyer  plate,  at  a  free  surface  of  the  sample. 

By  considering  the  physics  of  energy  deposition  and  shock  propagation  in  a 
porous  medium,  we  have  shown  that  conventional  techniques  do  not  yield  valid  measure¬ 
ments  of  the  true  solid-state  Gruneisen  parameter,  i.  e. ,  the  relation  between  pres¬ 
sure  and  specific  mergy  at  constant  volume,  because  on  a  microscopic  scale  de¬ 
termined  by  particle  size,  stress  relief  occurs  before  the  final  state  is  reached. 

This  is  true  for  the  case  of  shock  introduction  and  energy  deposition.  In  the  case  of 
energy  deposition,  simple  models  which  attempt  to  consider  the  result  of  stress 
relief  within  deposition  time  are  presented.  While  the  models  are  inadequate,  they 
serve  as  a  starting  point  for  understanding  the  physics  and  for  correlating  the  energy 
deposition  and  shock  introduction  cases. 

Energy  deposition  experiments  of  Shea,  Mazalle,  and  Avrami  (14),  and  shock 
introduction  experiments  of  Boade  are  analyzed.  In  addition,  we  have  applied 
some  elementary  considerations  of  a  porous  material  as  a  diffractive  medium  to 
comment  upon  the  difficulty  of  utilizing  macroscopic  experiments  for  investigating 
microscopic  parameters.  The  difficulty  being  that  local  unidirectional  one-dimensional 
strain  cannot  be  assumed  to  hold.  In  this  regard,  we  considered  some  flyer  plate 
work  in  polyurethane  as  reported  by  Butcher  (21). 

In  the  Introduction,  we  have  considered  the  importance  of  porous  materials 
from  an  academic  and  f  t>m  an  engineering  point  of  view.  Porous  materials  can  be 
used  for  shock  nritigatora.  and/or  as  thermal  standoff  materials.  Consequently,  the 
U.  S.  Army  and  the  shock  wave  community  would  like  to  understand  as  much  as  pos¬ 
sible  about  porous  solids.  Considerable  more  work  is  needed. 
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L  INTRODUCTION 


Porous  materials  are  any  materials  that  happen  to  have  a  mass  density  less 
than  the  maximum  possible  equilibrium  mass  density  at  the  pressure  and  tempera¬ 
ture  in  question.  Examples  of  such  materials  are  foams,  e.  g. .  polyurethane,  lees 
than  normal  density  explosives  and  porous  metals,  e.  g. ,  porous  tungsten  (2). 

Porous  materials  are  of  both  engineering  and  academic  Interest.  Polyurethane 
foams  are  utilized  as  shock  mitigators  by  having  the  shock  energy  absorbed  in  the 
process  of  compacting  the  foam  to  its  normal  solid  density.  Porous  metals  have 
engineering  applications  but  are  also  of  extreme  academic  importance.  If  a  solid 
of  normal  density  is  shock-loaded,  its  final  state  lies  in  a  single  curve  (the  Hugoniot 
curve)  in  the  PV  plane.  On  the  other  hand,  when  a  solid  such  as  porous  tungsten 
is  shock-loaded,  the  final  state  reached  will  reside  on  one  of  a  continuum  of  Hugoniot 
curves,  depending  upon  the  initial  porosity;  the  irreversible  processes  which  occur 
during  compacting  yield  PVT  states  for  the  compacted  solid  different  from  those  oc¬ 
curring  in  the  shock-loaded  material  having  initial  normal  density.  Thus,  by  shock¬ 
loading  an  initially  porous  material  to  bevond  complete  compaction,  it  is  possible 
to  investigate  broad  regions  in  PVT  ap? :  2.  This  academic*  property  of  porous  ma¬ 
terials  was  first  pointed  out  by  Zeldovich  (3). 

Another  important  application  of  porous  materials  is  in  the  area  of  explosives. 

It  is  well  known  that  the  initiation  sensitivity  of  explosives  to  mechanical  shock  depends 
upon  the  final  density  to  which  a  granular  explosive  may  be  pressed  ^).  This  area 
also  has  both  academic  and  engineering  applications. 

Soviet  and  Western  interest  in  inert  porous  materials,  judging  from  me  open 
literature,  has  been  in  different  pressure  regimes.  Western  interest +  (3, 6, 7, 8, 9)  has 
been  strong  in  the  area  of  shock  propagation  in  the  pressure  range  at  which  compaction 
occurs.  Or.  the  other  hand,  the  Russians  appear  to  have  their  main  effort  (2»  H)  in 
the  pressure  range  beyond  which  compaction  has  already  occurred.  Among  the  interesting 
ingredients  of  high-pressure  physics  which  the  Russians  nave  studied  via  porous  ma¬ 
terials  is  the  electronic  contribution  to  the  Griineisen  constant  (2). 

There  are  two  general  ways  by  which  porous  materials  may  be  examined  experi¬ 
mentally.  In  one  method,  energy  is  dumped  into  the  material  in  a  very  small  time, 
more  or  less  uniformly  as  a  function  of  position  and  the  shock  effects  accompanying  the 
resulting  thermal  expansion  may  then  be  studied.  A  typical  procedure  for  dumping  in 
the  energy  utilizes  an  electron  beam  machine  (12)  which  produces  (typically),  Mev 


+The  references  m  this  report,  for  the  most  part,  are  meant  to  be  illustrative  rather 
than  exhaustive. 
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electrons  in  large  numbers  within  a  pulse  width  from  20  nanoseconds  to  70  nanoseconds. 
The  elections  are  scattered  by  the  solid  parts  of  the  porous  material,  resulting  in  as 
increase  in  the  internal  energy  of  die  material.  The  second  general  method  involves 
the  impacting  of  the  porous  material  with  a  flyer  plate  (*3)  of  known  mechanical  prop¬ 
erties.  of  the  abode  originating  within  the  material  as  in  the  energy  dumping 

method,  the  shock  now  originates  at  die  interface  between  the  flyer  and  the  sample 
target. 


In  this  report,  we  wish  to  specifically  consider  an  aspect  of  some  electron  beam 
experiments  (M)  in  porous  samples  of  the  explosive  PETN  (Pentaerythritol  Tetra- 
nitrate).  Some  other  experimental  work  will  also  be  considered. 
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IL  THE  GRUNEISEN  CONSTANT 


There  are  two  possible  definitions  for  the  Gruneiaen  parameter,  one  micro* 
scopic,  and  the  other  macroscopic.  The  microscopic  involves  a  model  for  inter¬ 
actions  on  an  atomic  level,  while  the  macroscopic  definition  involves  measurable 
parameters. 

In  the  macroscopic  definition,  the  Griineisen  constant,  T,  is  defined  by: 


where  p  is  pressure,  p  is  mass  density,  and  E  is  internal  energy4.  Eq.  (1)  says 
that  if  one  deposits  thermal  energy  in  a  time  small  compared  to  the  time  necessary 
(defined  below)  for  significant  mass  motion,  then  the  pressure  In  the  solid  changes. 

The  quantity  T  determines  the  change  in  pressure,  and  is  a  property  of  the  structure 
of  the  material. 

The  time  criteria  mentioned  above  is  usually  taken  as  the  time  necessary  for 
an  acoustic  wave  to  traverse  an  electron  scattering  mean  free  path  in  the  solid.  Let 
such  a  time  be  tj.  If  the  pulse  duration  is  T,  then  the  constant  volume  definition  of 
Eq.  (1)  requires: 

T  >  ti.  (2) 

Physically,  this  means  if  Eq.  (2)  holds,  then  any  energy  inhomogeneitiea  introduced 
by  the  electron  beam  are  not  relieved,  by  acoustic  signals  before  the  energy  is  de¬ 
posited.  In  other  words,  the  energy  can  be  considered  as  deposited  instantaneously8 
if  Eq.  (2)  holds. 

4 Actually,  the  Griineisen  parameter  is  a  tensor  quantity,  and  the  stress  tensor  should 
be  substituted  for  the  pressure  in  Eq.  (1).  Such  an  anisotropy  effect  could  be  im¬ 
portant  in  this  report  where  we  are  interested  in  a  pressure  range  comparable  to 
the  yield  strength  of  the  solid  particles  of  PETN.  Perhaps  more  important,  the  pres¬ 
sure  range  of  interest  is  comparable  to  the  pressure  used  in  pressing  the  PETN  to 
an  initial  density.  Such  pressing  undoubtedly  introduces  anisotropy  into  the  material 
which  could  conceivably  show  up  in  the  experiment.  In  the  Shea,  Mazella,  and  A v rami 
experiment  l1*),  peak  pressures  generated  were  in  the  one-kilobar  range,  while  pres¬ 
sing  pressures  were  in  the  one  to  teu-kilobar  range.  Because  of  a  lack  of  sufficient 
experimental  evidence  to  decide  this  question  at  the  present  time,  we  will  assume 
here  that  Eq.  (1)  is  correct  in  its  present  isotropic  form. 

For  those  computer  codes  which  can  handle  energy  deposition  and  stress  propaga¬ 
tion  simultaneously,  T  should  be  replaced  in  Eq.  (2)  by  the  time  cycle.  If  the  energy 
deposition  is  by  other  than  electrons,  the  appropriate  mean  free  path  must  be  used. 
Since  this  report  makes  a  point  of  physics,  rather  than  technique,  the  possibility  of 
handling  hydrodynamics  and  energy  deposition  simultaneously  is  delayed  until  the 
discussion  section  of  the  report. 


One  o*  the  main  points  of  this  report  is  that  the  time  criteria  of  Eq.  (2) 
is,  in  general,  not  valid  for  a  porous  material.  Ai!  it  takes  to  invalidate  that  cri¬ 
teria  is  that  Ae  average  dimension  of  the  solid  particles  making  up  the  material  be 
small  compared  to  the  appropriate  radiation  mean  free  path.  For  such  small  par¬ 
ticles,  the  eaergy  deposition-induced  thermal  stresses  are  relieved  by  thermal  ex¬ 
pansion  prior  to  completion  of  the  energy  deposition  process.  Table  I  below  lists 
some  electron  energy  absorption  mean  free  paths  for  a  variety  of  materials  and 
energies. 


TABLE  I 

ELECTRON  MEAN  FREE  PATH  AS  A  FUNCTION  OF  ENERGY 

Material 

Mean  Free  Path 

Electron  Energy 

Al 

100  n 

0.4  Mev 

PETN 

500  fi 

0.4  Mev 

Be 

1,000  fi 

0. 6  Mev 

Quartz  Crystal 

4, 000  n 

2  Mev 

PETN 

1,000  ft 

4  Mev 

Lead  Styphnate 

1,000  n 

4  Mev 

As  contrast  to  the  typical  dimension  of  10° n ,  i.  e. ,  0. 1  cm,  of  Table  I,  an 
acoustic  wave  propagating  at  3  x  10^  cm/sec  can  relieve  the  stress  to  the  center  of 
a  18 Op  particle  in  30  nanoseconds.  Consequently,  an  electron  beam  machine  having 
a  pulse  width  of  30  nanoseconds  (FWHM)  would  not  be  able  to  utilize  the  criteria  of 
Eq.  (2)  for  experiments  with  porous  solids  having  particle  sizes  of  18Cp  or  less.  It 
is  interesting  to  note  that  the  30-nanosecond  figure  is  typical  of  the  faster  machines 
in  use  today.  In  the  work  by  Shea,  Mazella,  and  Avrami  (hereafter  called  SMA),  the 
average  grain  size  of  the  PETN  was  observed  to  be  250u>  The  electron  beam  was  at 
0. 4  Mev  with  a  pulse  width  of  25  nanoseconds  (FWHM),  Thus,  if  the  PETN  in  question 
was  porous  enough  so  that  a  significant  number  of  grains  were  not  in  intimate  contact 
along  a  grain  boundary,  then  the  time  criteria  of  Eq.  (2,  would  be  invalid. 

We  hasten  to  say  that  the  work  of  SMA  is  still  interesting  in  that  it  yields 
the  magnitude  of  the  pressure  pulse  resulting  from  a  thermal  energy  deposition.  In¬ 
deed,  the  experimenters  were  aware  of  the  point  with  respect  to  time  criteria  made 
above  at  the  time  of  their  experiment,  and  that  is  why  they  speak  in  terms  of  an  ''ef¬ 
fective"  Griineisen  constant  in  their  paper.  The  object  now,  however,  is  to  see  whether 
the  measurements  made  can  yield  any  useful  information  concerning  microscopic 
processes.  Table  II  lists  the  pertinent  results  of  the  SMA  experiments. 
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TABLE  na 


-  - - — -  —  — 

ENERGY  DEPOSITION  PROPERTIES  OF  PETN 

b 

Poo  4>0 

Density 

gm/cm** 

Sound  Speeds  g 
cm/sec  tc  10**° 

Effective 

r 

0.  95 

1. 67 

2.8 

1.2 

1.59 

2. 4  ±  0. 1 

0.51 

0.87 

1. 54 

1.  8  ±  0. 3 

0.15 

0.84 

1. 48 

1. 7  ±  0. 3 

(0.  07  to  0. 23) 

aThis  table  was  taken,  with  minor  changes,  from  reference  14. 
b  pQ  denotes  the  normal  density.  poa  initial  porous  material  density. 


Although  SMA  report  an  experimental  error  in  the  measured  T  of  20  percent 
or  more,  the  trend  as  a  function  of  density  is  clearly  discernible. 
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in.  MICROSCOPIC  EFFECTS  AND  THE  GRUNELIEN  CONSTANT 


In  this  section,  we  discuss  some  of  the  physics  which  might  be  important 
in  interpreting  the  measured  or  effective  Griineisen  constant  Integrating  Eq.  (1) 
gives* 


P=pir(E-Ei)  (3) 

where  pi  denotes  the  initial  density  at  which  energy  is  deposited  in  the  constant  volume 
approximation,  T  is  assumed  independent  of  energy  density,  and  E}  is  the  specific 
energy  density  prior  to  energy  deposition.  That  T  is  indeed  independent  of  E  has 
been  borne  out  experimentally  by  the  success  of  the  linear  P  vs  E  curve  in  fitting 
the  data  (*4). 


There  are  two  relevant  microscopic  effects  associated  with  the  energy  deposi¬ 
tion  aspect  of  porous  solids.  First,  there  is  the  possibility  that  surface  energy  ef¬ 
fects,  because  of  the  large  internal  surface  area  expected  in  a  porous  solid,  can  play 
a  role  in  determining  the  pressure,  P.  Second,  is  the  expected  and  mentioned  stress 
relief  in  the  solid  particles  of  the  porous  material. 


Let  us  write  for  the  surface  energy,  yfl: 
apY 


=  20 


(4) 


where  aQ  is  the  lattice  parameter,  and  Y  is  Young's  modulus  for  the  solid  particles. 
Using  values  for  a  typical  solid: 


Y  =  1012  dynes/cm2 ,  Sq  =  2  x  10“8  cm, 

we  find  that: 

v  =  103  ergs/cm2. 

®  i 

The  heat  of  formation,  AHf  =  5  x  10’  cal/gram  formula  weight.  For  a  solid  of  pQ  - 
5  gm/cm3  and  a  gram  formula  weight  of  50,  we  have  an  E0  of: 

E0  =  5  x  103  cal  ~  2  x  1011  ergs 


+In  Eq.  (3),  we  assume  that  the  thermal  pressures  generated  are  considerably  larger 
than  the  pressure  before  deposition,  i.  e. ,  atmospheric  pressure,  and  also  larger 
than  anyintemal  stresses  generated  in  the  process  of  pressing  the  porous  material 
to  its  initial  density.  There  is  a  possibility  that  this  last  assumption  is  invalid 
since  typical  measured  pressures  (SMA)  were  in  the  fraction  of  a  kilobar  range, 
while  pressing  pressures  were  in  the  one  to  ten-kilobar  range. 
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for  one  cm3  of  a  typical  solid.  If  the  same  cm^  of  material  is  now  considered  to 
have  103  cm2  of  internal  surface  area,  S,  (a  factor  of  103  more  chan  external  surface 
area): 

y8S  =  106  ergs. 

Comparison  of  E0  and  ysS  would  seem  to  indicate  that  surface  energy  effects  can 
be  neglected.  This  view  is  further  supported  by  the  idea  that  surface  energy  is  not 
expected  to  change  during  energv  deposition;  thermal  expansion  should  allow  stresses 
built  up  by  deposition  to  be  more  or  less  instantly  relieved  so  that  neglecting  tem¬ 
perature,  the  surface  remains  in  a  constant  state. 

Let  P}  denote  the  pressure  in  the  normal  density  solid,  and  P2  the  pressure 
in  the  corresponding  porous  solid  of  initial  density  p0o: 


P1  s  Po  *1  (E  -  Eq) 

(5a) 

P2  5  Poo  r2  <E  ‘  Eo> 

(5b) 

Since  it  has  been  shown  that  surface  energy  considerations  are  not  important,  the 
same  specific  internal  energy  has  been  used  in  Eqs.  (5a)  and  (5b).  Now  let  Vs  be 
the  volume  occupied  by  the  solid  particles,  and  Vp  that  occupied  by  the  voids.  We 
wish  to  investigate  the  relation  between  and  ^  as  a  function  of  various  models 
for  averaging  the  pressure  over  the  solid  particles  of  the  porous  material.  The 
averaging  takes  place  with  the  mass  motion  frozen  at  the  time  at  which  energy  deposi¬ 
tion  ends.  One  of  the  simplest  averaging  models  is  given  by: 


(6) 


where  simply  denotes  the  pressure  with  this  particular  averaging  model.  Eq. 

(6)  says  that  the  pressure  is  reduced  from  that  of  the  normal  density  solid  by  being 
smeared  out  over  the  entire  volume  to  the  solid  particle.  If  m  is  the  mass  of  the  solid 
particle: 


Although  the  model  given  by  Eq.  (6)  is  unphysical,  it  is  still  interesting  to  note  that 
it  takes  an  exrx)nent  n  greater  than  ten  to  fit  the  data  of  Table  n,  while  n  =  1  yields  a 
Griinoisen  parameter  independent  of  porosity. 
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CO  |£>- 


The  next  averaging  model  we  consider  is  more  physical  and  is  shown  in  fig¬ 
ure  1  following: 


FIGURE  1:  An  Averaging  Model  for  Effective  Pressure  Due  to 
Energy  Deposition  in  a  Porous  Solid 


Averaging  the  pressure  over  the  solid  particles: 


where  we  are  coniiderittg  the  particles  to  be  spherical  in  Eq.  (10).  Using  - 
tr  as3  =  m/po,  we  thus  find: 


+V 


(V  +  V  ) 
'  s  p' 


(11) 

(12) 

(13) 
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The  model  leading  to  Eq.  (13)  is  equivalent  to  assuming  that  a  relief  wave  propagates 
toward  the  center  ol'  the  particle,  and  that  the  wave  decreases  in  amplitude  as  it  propa¬ 
gates  toward  the  particle  center.  It  would  also  appear  that  a  further  assumption,  that 
the  wave  just  reaches  the  center  of  the  particle  at  the  end  of  the  energy  deposition,  has 
been  made.  Once  again,  the  averaging  model  is  unphysical  but  interesting.  We  find 
an  effective  Griineisen  parameter  which  is  dependent  upon  the  spherical  geometry,  but 
independent  of  porosity.  If  anything,  we  have  shown  that  the  naive  approaches  will 
not  work. 

We  feel  that  this  type  of  work  should  be  carried  further.  A  more  powerful  and 
more  physical  approach  than  that  presented  above  would  use  a  computer  program  for 
stress  wave  analysis.  The  program  would  have  to  handle  energy  deposition  and  hydro¬ 
dynamics  at  the  same  time  i^).  Having  such  a  code,  the  proposal  is  then  to  look  at 
the  pressure  time  distribution  in  a  slab  of  solid  material,  i .  e. ,  pQ  =  p  ,  where  the 
slab  is  so  thin  that  the  acoustic  transit  time  through  it  is  small  compared  to  the  time 
associated  with  energy  deposition. 

There  is  a  definite  lack  of  data  specifically  relating  to  Griineisen  constant  of  - 
porous  materials  before  compaction  occurs.  It  is  possible,  however,  to  analyze  ex¬ 
isting  shock  velocity  versus  particle  velocity  data  to  obtain  such  information.  Figure 
2  following,  is  based  upon  data  taken  from  Boade  (1S)  on  sintered  porous  copper.  If 
we  let  Us  be  shock  velocity,  and  up  be  particle  velocity,  then  the  extrapolations  (dashed 
lines)  in  figure  2  can  be  represented  by: 

Us  =  A  +  B  up  (14) 

Let  ro  be  Griineisen  parameter  corresponding  to  the  extrapolated  zero  pressure 
state.  It  is  easily  shown  0  9)  that: 

T0  =  2B  - 1  (15) 

It  is  then  a  simple  matter  to  apply  Eq,  (15)  to  figure  2. 

For  the  data  shown  in  figure  2,  experimental  points  only  exist  for  the  67. 7 
percent,  82.  8  percent,  and  solid  density  sample.  The  other  curves  are  theoretical. 
Applying  Eq.  (15)  to  the  mentioned  data  gives  the  results  shown  in  Table  III. 
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Shock  Velocity,  mm/p  sec 


Solid  Copper 

p  Q  =  8. 93  g/cm3 

U  =  4. 02  +  1. 46  u_ 
8  F 


/  / 
/  99% 


95%  * 

'/ 

/ 


(/  82.9%  '' 

'  A'V 

67.7%^  - 


Particle  Velocity,  mm/p  sec 

FIGURE  2:.  Shock  velocity  vs  particle  velocity  for  sintered  porous 

copper  (from  Boade13).  Percentage  porosity  is  indicated 
by  Poo ^  anc*  dashed  lines  are  meant  to  extrapolate  to 
zero  pressure  data. 
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!  TABLE  IH 

j 

i  Criioeisea  Parameter  for  Sintered  Porous  Copper 

i 

poc/(>c 

0.677 

2.79 

0.829 

4.67 

1.00 

2.25 

The  data  thus  show  that  the  Grlineisen  parameter  decreases  as  the  degree  of 
porosity  increases*  This  is  similar  to  the  averaging  model  of  Eq.  (8).  On  the  other 
hand,  the  F0  value  corresponding  to  the  solid  density  copper  is  lower  than  that  of 
aay  of  the  experimentally  observed  porous  density  F0  values.  This  last  point  is  both 
unexpected  and  interesting.  It  correlates,  however,  with  an  observation  made  by 
SMA  (14)  that  the  calculated  Grineisen  parameter  for  normal  density,  i.  e. ,  solid, 
PETN  is  0.  8  (with  perhaps  20  percent  error),  which  is  smaller  than  theT  associated 
with  Poo/po  ~  0. 95,  as  shown  in  Table  IL  For  the  moment,  this  phenomenon  is  not 
understood. 
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IV.  SHOCK  PROPAGATION  IN  POROUS  MATERIALS 


Id  this  section  we  would  like  to  discuss  the  relationship  between  the  Griineisen 
data  obtained  from  the  impact  experiments  of  Boade,  and  the  energy  deposition  data 
obtained  by  SMA.  We  will  then  further  discuss  an  aspect  of  the  physics  associated 
with  shock  propagation  in  porous  materials. 

Consider  a  porous  material  modeled  as  shown  in  figure  3: 


Shock 

- 

Propagation 

Direction 


FIGURE  3:  Model  for  a  Porous  Solid.  (The  shaded  areas  represent 
solid  particles,  and  the  unshaded  areas  are  voids. ) 

If  a  shock  of  amplitude  too  small  for  compaction  propagates  into  the  material  shown 
in  figure  3?  we  first  expect  something  akin  to  the  Hugoniot  state  of  a  normal  density 
solid  to  be  transmitted  across  the  contacting  surfaces  of  the  solid  particles.  The 
pressure  corresponding  to  such  a  Hugoniot  state  would  then  be  relieved  in  time  by 
rarefactions  entering  the  solid  particles  through  their  free  surfaces.  For  solid  par¬ 
ticles  with  diameter  small  compared  to  the  width  of  the  shock  front,  the  final  state 
would  thus  be  a  function  of  the  Hugoniot  state  in  the  normal  state  plus  a  relief  factor. 


Since  the  Us  versus  up  relation  is  measured  for  the  final  state,  and  thus  so 
is  the  Griineisen  relation,  the  degree  of  porosity  must'  play  a  role  in  determining  F. 
Greater  porosity  means  greater  relief  of  the  normal  density  pressure  state,  a  re¬ 
duced  5US  value,  and  thus  a  reduced  T.  The  same  physical  phenomenon  appears 


to  be  occurring  in  the  shock  propagation  determination  of  F  as  occurs  in  the  energy 
deposition  case. 
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MICROSCOPIC  EFFECTS  ON  SHOCK  PROPAGATION 


V. 


An  observer  riding  on  a  shock  front  in  a  porous  material  such  as  shown  in 
figure  3  would  see  a  succession  of  diffraction-grating  planes.  The  effect  of  diffraction 
in  general  is  to  change  the  geometry  of  the  wave  front.  In  our  case,  the  change  is 
from  plane  wave  geometry  to  nonplane  wave  geometry.  Further  specializing,  we  can 
say  that  on  a  microscopic  particle  size  level,  the  mass  disturbance  is  not  unidirection 
one -dimensional  strain*  even  though  it  may  be  unidirection  one-dimensional  strain  on 
3  macro4. co  pic  level.  One-dimensional  strain  occurs  when  the  dimensions  of  the 
sample  perpendicular  to  the  direction  of  shock  propagation  are  large  compared  to  the 
sample  dimensions  parallel  to  die  direction  of  shock  propagation.  If  the  relative  size 
oi  the  dimensions  are  reversed,  then  the  propagation  is  classified  as  one-dimensional 
stress.  On  a  microscopic  level,  for  spherical  particles  with  some  void  separation 
between  them  as  shown  in  figure  3,  we  thus  expect  a  change  in  the  geometry  of  the  wave 
front. 


If  one  is  only  interested  in  characte.iring  the  macroscopic  shock  properties  of 
a  porous  solid,  then  the  unidirection  one-d  unentional  strain  view  of  the  problem  (with 
appropriate  perpendicular  and  parallel  dimensions)  is  valid.  If,  however,  one  wishes 
to  use  the  data  from  the  macroscopic  experiment  to  determine  microscopic  parameters 
(such  as  a  relaxation  time),  large  errors  could  result.  Physically,  the  reason  for 
such  errors  is  that  diffraction  modifies  the  way  in  which  the  principle  of  the  conserva¬ 
tion  of  momentum  is  applied. 

In  figure  4  following,  we  illustrate  the  result  of  microscopic  diffraction.  Or¬ 
iginally,  a  plane  wave  is  propagating  ir.  the  z  direction.  After  diffraction  has  occurred, 
small  segment  of  the  wave  front,  with  area  (#3)*  moves  in  the  £  direction  (normal  to 
the  £  axis).  s 


*By  one-.iimcTisnuial  strain  we  imply  plane  wave  propagation  wiili  particle  motion  only  in 
the  direction  of  the  wave  propagation.  The  one-dimensional  strain  case  is  the  typical 
configuration  in  shock  wave  theory  and  experiment. 
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For  brevity,  assume  that  the  local  particle  motion  is  also  in  the  4  direction. 

This  simplification  neglects  the  possibility  of  local  transverse  waves  being  generated 
due  to  the  porous  structure  of  the  medium,  ii  other  words,  we  have  local  one- 
dimensional  strain  in  other  than  the  z  direction.  If  rr>  represents  the  momentum  flux  (2®) 
in  the  4  direction,  then  (2®) :  * 

(AS)^  s(P+pu^u^)  (AS)^  ,  (16) 

where  p  Is  pressure,  p  is  mass  density,  and  u*  the  particle  velocity  in  the  4  di¬ 
rection. 


CAS)^  =  (AS)Z  cos  8  (17) 

Let  Pz  be  the  momentum  transferred  per  time  across  (A  S)2. 

•  • 

Pz  3  nz  (AS)Z  =  P^  cos  0  (18) 

Pz  =  (AS)^  cos  0  (19) 

trz  =  ( P  +  p  u^  u^  )  cos2  0 ,  (20) 

or 

nz  -  p  cos2  0  +  p  uz  uz.  (21) 

Two  points  now  need  to  be  made.  Fi  rst,  when  employing  the  conservation 
of  momentum  in  shock  wave  problems,  people  usually  equate  the  right  hand  side  of 
Eq.  (21)  on  each  side  of  the  shock  front  (with  0=0).  Second,  even  though  no  net  de¬ 
viation  from  macroscopic  one-dimensional  strain  in  the  z  direction  is  observed,  a 
nonzero  average  value  for  cos2  0  is  expected  to  exist. 
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Butcher’21'  uses  Eq.  (20)  with  0  =  0  tc  investigate  abode  propagation  in 
0. 96  gm/'cc  poiyurethane  foam,  he  was  interested  in  the  relaxation  time,  Tq,  for  void 
collapse  (closure;  and  generated  an  equation  of  die  form: 


„=A(g-Ctq)/T0 
(B  +  Co) 


(22a) 


where  A,  B.  and  C  may  be  regarded  here  as  constants.  Oir  point  is  that  an  equation 
more  in  the  form  of:  _ 

A  (a  -  cr„)  cos20/to 

u2  =  V  (22b) 

(B  +  Cff  cos20) 

should  have  been  used.  In  Eqs.  (22),  or  denotes  stress,  or^  an  equilibrium  stress 
value,  and  cos20  an  average  value  of  cos20.  Such  a  correction,  at  least  in  the  nu¬ 
merator  of  Eq.  (22b),  would  result  in  a  larger  value  of  r0  than  previously  estimated. 

A  larger  value  of  r0  would  be  in  keeping  with  some  other  ideas  (based  upon  a  "sim¬ 
plified  gap  model")  put  forward  by  Butcher  At  the  moment,  we  do  not  have  a 
good  model  for  calculating  cos-0. 
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